A precise CNOT gate in the presence of large fabrication induced variations of the 

exchange interaction strength. 
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We demonstrate how using two-qubit composite rotations a high fidelity controlled-NOT (CNOT) 
gate can be constructed, even when the strength of the interaction between qubits is not accurately 
known. We focus on the exchange interaction oscillation in silicon based solid-state architectures 
with a Heisenberg Hamiltonian. This method easily applies to a general two-qubit Hamiltonian. We 
show how the robust CNOT gate can achieve a very high fidelity when a single application of the 
composite rotations is combined with a modest level of Hamiltonian characterisation. Operating the 
robust CNOT gate in a suitably characterised system means concatenation of the composite pulse 
is unnecessary, hence reducing operation time, and ensuring the gate operates below the threshold 
required for fault-tolerant quantum computation. 

PACS numbers: 03.67.Lx, 82.56.Jn, 85.30.De 
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Introduction 



The ability to correct errors arising from the construc- 
tion or operation of any quantum computing architec- 
ture is essential for a successful implementation. With- 
out the ability to correct the random and/or systematic 
errors that arise throughout operation, the implemen- 
tation of large scale quantum algorithms is hopelessly 
undermined. In a realistic device the threshold for fault- 
tolerant quantum computation is likely to be well below 
10~ 4 , placing severe constraints on the tolerable mag- 
nitude of errors due to decoherence or lack of precision 
in quantum control. This work focuses on minimising a 
particular type of systematic error, namely, uncertainty 
in the coupling strength of two-qubit devices as a result 
of imperfect fabrication, which causes systematic under- 
or over-rotations. We use recently developed two-qubit 
composite rotations to correct for this uncertainty in the 
strength of the electron spin exchange interaction in Si:P 
based architectures^ 2 -. Our results also apply more gen- 
erally and could be used to correct this type of systematic 
error in a range of solid-state systems. 

The strength of the exchange interaction coupling be- 
tween donors in silicon based solid-state architectures is 
known to be highly sensitive to donor placement. The 
cause of this is the inter-valley interference between the 
six degenerate conduction band minima of silicon, result- 
ing in oscillations of the exchange coupling strcngth2ii£i&. 
Exact positioning of donors to better than 2-3 sites is 
difficult^ and therefore we expect significant uncertainty 
in the un-biased strength of the coupling between donors. 
The uncertainty in our knowledge of the coupling, leads 
to error in gate operation. Systematic errors of this 
kind are correctable using composite rotations. Experi- 
mental applications already exist in a variety of quan- 
tum systems demonstrating the usefulness of compos- 
ite rotations for ensuring robust operation s 8 i 9 i 10 ' n i 12 ' 13 . 



Recently two-qubit composite rotations have been con- 
sidered for systems with uncertainty in their coupling 
strength!^. 



In this paper, we follow the method for creating a ro- 
bust controlled-NOT (CNOT) gate developed in Ref. GJ 
and quantitatively study the performance of the robust 
CNOT gate using simulated exchange oscillation data. 
We specifically consider the global Si:P electron spin con- 
trol case where the interaction is of Heisenberg type and 
gate times are in the 0(10-100 ns) regime. This tech- 
nique is readily generalisable to any two-qubit Hamilto- 
nian, and for a full treatment, the reader is directed to 
Ref. M 



Misplacement of donors by only one implantation site 
can lead to large variations in the exchange coupling 
strength, even in Si:P systems with voltage bias applied 
to top gates^, meaning a single application of the com- 
posite rotations may not be enough to guarantee a high 
fidelity CNOT gate. Concatenating the pulse by feed- 
ing it back into itself can help to achieve correction to 
a higher level, however, performing multiple concatena- 
tions costs a large increase in time. In certain cases using 
composite rotations alone will not improve the fidelity of 
the operation above an uncorrected CNOT gate, as the 
composite rotations are designed to work within a specific 
uncertainty range. Wc show that in unison with Hamilto- 
nian characterisatio n 16 ! 1 7 , the process of experimentally 
determining a Hamiltonian, a single application of the 
composite rotations guarantees a high fidelity CNOT op- 
eration with an error rate below the fault-tolerant error 
threshold. Operating the CNOT gate this way helps re- 
move the need for concatenation, and strikes a balance 
between fully characterising the system and using com- 
posite rotations to construct robust operations. 
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FIG. 1: Circuit diagram for a CNOT gate constructed from 
an Ising interaction, where H is a Hadamard gate, Z_n. = 
exp [i\az) and Ui (f ) = exp {-i\<j z ®oz)- 



I. CONSTRUCTING ROBUST GATES USING 
COMPOSITE ROTATIONS 

Composite rotations have been widely used in NMR 
experiments to correct for pulse length errors and off- 
resonance effect o 18 ' 19 . In the case of pulse length errors, 
a deviation of the field strength from its nominal value 
leads to systematic under- or over- rot at ions. Although 
originally designed for applications involving single spin 
quantum systems, composite rotations may be extended 
to two-spin operations. In the context of quantum com- 
putation, only a certain class of composite rotations, 
sometimes referred to as fully compensating pulses are 
applicable, as they work on any initial state. Using these 
fully compensating pulses, the application of composite 
rotations for constructing robust two-qubit gates against 
pulse length error has already been found for an Ising 
Hamiltonian 1 ^, and a general two-qubit Hamiltonian 1 ^. 

In Ref. [TH, it was noted that for a general two-qubit 
Hamiltonian expanded in the Pauli basis, 



H 



E 



(1) 



any interaction term can be effectively extracted using 
a technique called term isolation^. The isolation of a 
given term will in general not be exact but can be made 
arbitrarily accurate. This result is particularly useful and 
can be used to isolate the Ising coupling term, Jzz, such 
that we can construct a CNOT gate from this interaction 
as in Fig.[TJ In the case of the Heisenberg interaction with 
isotropic couplings, 

Hn = J(a x <%> ax + ay ® a Y + az ® a z ), (2) 
the isolation of the Jzz term is exact, 



exp (-iJzztvz <8 a z ) 



---{Z v ® I) exp (-iHnt) 
x (Z v <E> I) exp (-iH H t) , (3) 



where for single qubit gates Z a is a rotation about the 
oz axis by an angle a, and similarly for other operators, 
Jzz = 2 J, and the global phase factor is included. 

We now consider constructing a robust CNOT gate 
using composite rotations, whereby we replace the inter- 
action term with one created using composite rotations. 



Doing this compensates for any uncertainty in our knowl- 
edge of the exchange interaction coupling strength, J. In 
Fig. [5] the entire process of constructing a robust CNOT 
gate from composite rotations is demonstrated schemat- 
ically. 

In an ideal system with a perfectly characterised cou- 
pling strength, the evolution operator generated by the 
Ising interaction is 



Oo = m (0) = exp 





' l 2 az ® ° z 



(4) 



Here, 0q is a two-qubit rotation by an angle about the 
az ® oz axis. In general, a is a two-qubit rotation by 
an angle around an axis tilted from the az ® az axis 
towards the az ® ax axis by an angle a, 



exp 



az cos a + az ® ax sin a) 



This two-qubit rotation is achievable via, 

0a = (i®Y a )9 o (i®y_„). 



(5) 



(6) 



We make the assumption that all single qubit unitaries 
are error free, but note that single qubit operations may 
also be made robust using existing techniques developed 
in the context of NMR. 

In reality a fractional error, A, in the two-qubit opera- 
tion will be present due to the uncertainty in our knowl- 
edge of the actual coupling strength, Jzz, 



A = 



Jzz 

~jV 



- 1. 



(7) 



Here, Jp is our prediction of the Ising coupling strength 
based on the targeted donor separation. Therefore the 
actual rotation performed will be 



i 0) = u(e) 



exp 



The superscript of 0a in the above equation indicates 
the implementation level of the actual (non-ideal) rota- 
tion, with "(0)" being an uncorrected implementation 
and higher levels signifying subsequent corrections from 
composite rotations. The implementation level should 
not be confused with concatenation level, (e.g., 2 nd im- 
plementation level is the 1 st concatenation level). 

It has been previously noted that single qubit com- 
posite rotations can be extended to two-qubit compos- 
ite rotations for use in quantum computatio n 14 ' 15 using 
fully compensating pulses. A class of these composite 
rotations known as BB l 19 ' 21 is particularly useful for ap- 
plications involving quantum computation^ 2 -. Replacing 
the pulse O °^ with the symmetrised BB1 class composite 
pulse 



-i-(l + £)a z ®a z 



(8) 



9 (i) 



{0/2)f 



4 0) <M 0) 



(0) 



(9) 
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WITHOUT CONCATENATION 



Start with 

V , J 



Isolate Ising component 

exp (-iJzztvz ® <rz) = -(ZV ® I) exp {-iH H t) (Z w ® /) exp (-iH K t) 
V , J 



Construct composite rotations 

e a) = W2) (o) ir (o) 27r (o) 7r (o, {0/S) m 



Construct robust CNOT gate 
CNOT=(J ® H)(Z_ | ® Z_ | )$ <*>(/ ® H) 
for 2 = 1,2,... 



WITH CONCATENATION 




Correct <jz ® oz operations to higher level 
for i = 1,2, . . . 



Re-isolate Ising component 

e«* = [pr„ 55 x w ) (i ® z„) [e/(2iv r )]<* ) (/ ® z,,) [e/(2/v r )]^] 

for » = 1,2, ... 

where JV r is the level of re-isolation 



FIG. 2: Procedural flowchart for constructing a robust CNOT gate using composite rotations, and concatenating to higher 
implementation levels. 



where 4> = arccos(— O/Att), will result in a higher fidelity 
operation. The fidelity of an operation is defined as 



Tr 


~W (9) Ui (6)] | 


Tr 


U\ (9) U, (9) 



(10) 



We may re-isolate the Ising component Jzz again to ar- 
bitrary accuracy as in Fig. [2] The re-isolated Ising com- 
ponent can then be used to correct to even higher order 
by passing this pulse back into each of the constituents 
of Eq. [51 (see Fig. ^ . In principle there is no limit to 
how often this concatenation can be done, however, the 
increase in gate time means that in practice this process 
will be limited by the decoherence time of the system in 
which the CNOT gate is being implemented. In Fig.[3]the 
performance of the uncorrected CNOT gate is compared 
to the robust gate for various implementation levels, as 
originally calculated in Ref. [IH Notice each subsequent 
implementation level performs better over a larger range 
of the fractional error, A. 

We now apply the robust CNOT gate to the Si:P ar- 
chitecture with large fabrication induced variations (and 
hence uncertainty) in the exchange interaction strength. 



II. CORRECTING FOR AN UNKNOWN 
EXCHANGE INTERACTION STRENGTH 

Systematic errors arising from imperfections in the 
fabrication process arc correctable. In Kane type 




-1 -0.5 0.5 1 



A (fractional error) 

FIG. 3: CNOT error, (1 - T), as a function of the frac- 
tional error in our knowledge of the coupling strength, A, 
for various implementation levels. These composite rotations 
provide improvement over an uncorrected implementation for 
A G (-1,1). The fidelity of a CNOT gate constructed from 
the Heisenberg interaction using composite rotations was orig- 
inally calculated in Ref. Il5l . 



architectures^ where phosphorus donors are implanted 
into an isotopically pure 28 Si matrix, two fabrication pro- 
cesses are being pursued concurrently^. The top down 
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FIG. 4: Exchange couplings for donors at fee lattice sites mis- 
placed by a distance <5 in all directions from the target sep- 
aration of 20.634 nm (in the [100] direction). The exchange 
coupling strengths are given as a fraction of the target cou- 
pling strength, Jo. 



approach uses ion beam implantation of phosphorus ions 
incident on the silicon substrate. Precise placement of 
phosphorus donors is limited in this approach due to 
scattering off the silicon atoms, in a process known as 
straggling. State of the art top down fabrication re- 
sults in placement uncertainties of O(10 nm)2i. The 
bottom up approach offers atomically precise fabrication 
using a phosphine gas. The gas is applied to a hydro- 
gen terminated silicon substrate, where scanning tunnel- 
ing microscopy has removed individual hydrogen atoms 
from the hydrogen mono-layer at the desired implanta- 
tion sites. Once the phosphorus is integrated into the 
substrate, the mono-layer is removed and overgrown with 
silicon. Small deviations from target implantation by of 
0(1 nm) (approx. 2-3 sites) can still occur during the 
annealing process^. 

The exchange coupling J of the Heisenberg Hamilto- 
nian (see Eq.[2|), is highly sensitive to donor electron wave 
function overlap. This means that even small deviations 
from the targeted implantation sites can lead to large 
variations in the exchange coupling between donors^. 
Calculated variations in the strength of J for small de- 
viations from the targeted donor separation are shown 
in Fig. 2] This calculation was performed using the 
Heitler-London formalism, where the wave functions for 
the phosphorus donors in silicon were expressed in Kohn- 
Luttinger effective mass form, with Bloch states explic- 
itly computed using the pseudopotential fit to the band 
structure. Details can be found in Ref. 0- Importantly, 
this type of systematic error is correctable using the com- 
posite rotations described above. 

We consider donors to be implanted along the [100] 



FIG. 5: CNOT fidelity as a function of donor separation in 
the [100] direction for various implementation levels. The 
resulting fidelities are determined based on a target donor 
separation of 20.634 nm. Note that interpolating curves be- 
tween lattice sites indicate donor separation scenarios for a 
given implementation, and vertical dotted lines guide the eye 
between implementations. 

lattice direction, as oscillations are suppressed in this di- 
rection and the coupling is stronger relative to the [110] 
and [111] directions, meaning this is the preferred direc- 
tion for device fabrication. In an uncharacterised system 
we assume that the exchange interaction strength is Jo 
and will be determined by the target donor separation 
and bias on the control gates. Fabrication induced donor 
misplacement will cause the true exchange interaction 
strength, J, to be quite different from Jo. The fractional 
error in our knowledge of the coupling strength is 

A = ^-l. (11) 
Jo 

These composite rotations will only provide an improve- 
ment over an uncorrected implementation for |Arj| < 1. 
For |Ao| > 1 these composite rotations are actually 
outperformed by the uncorrected implementation, so if 
J > 2 Jo then the composite rotations provide a less ro- 
bust operation. Interestingly, provided J 7^ 0, in which 
case we have no entangling operation, composite rota- 
tions will correct for any J G (0, 2J ). 

Implementing the gate based on the target coupling 
strength Jo, the fidelity of the resulting CNOT opera- 
tion will be determined by the size of the fractional error 
Ao in the actual coupling strength. In Fig. Owe demon- 
strate the resulting CNOT fidelity for a number of donor 
separations in the [100] direction when the target separa- 
tion is 20.634 nm. The results show that using composite 
rotations improves the fidelity of operation for the CNOT 
gate. For example, if the actual separation is 21.720 nm, 
one application of the composite pulsing scheme improves 
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the fidelity from ~ 0.93 to ~ 0.99, whilst a second ap- 
plication brings the fidelity above 0.9999. The successive 
improvements due to the various levels of pulse concate- 
nation do however come at the expense of operation time. 
We examine this issue in the following sections. 

A. Gate count 

The robust CNOT gate outperforms the uncorrected 
CNOT gate given an error in the targeted coupling 
strength, J , for |A | < 1. Each level of concatenation 
provides further improvement, however the cost of this 
improvement is an exponential increase in the total num- 
ber of gates required. An unavoidable consequence of this 
is an increase in the time required to perform these ro- 
bust operations. To be of use for quantum computation 
we need to be able to perform many precise operations 
within the decoherence time of the system. Minimising 
the time taken to perform a robust CNOT gate is essen- 
tial. We show how this can be achieved in Sec. IIIII util- 
ising Hamiltonian characterisation. Below, we consider 
the actual time costs of concatenated composite pulse 
correction. 

An uncorrected CNOT gate requires only 6 single qubit 
gates and 2 two-qubit gates. In comparison, a raw gate 
count for the number of single qubit gates required in 
constructing the robust Ising interaction for the CNOT 
gate yields 

n\ = 16, 

m = HW r (rai_i+2) + 6, i = 2,3,..., (12) 

where m is the number of single qubit gates required 
for the i th implementation level, and N T , which we as- 
sume to be constant, quantifies how much we re-isolate 
the Ising term for pulse concatenation. Constructing a 
robust CNOT gate requires an additional 4 single qubit 
gates, such that the total number of single qubit gates 
required, n lq , is 

n] q = m + 4, * = 1,2, ... . (13) 

The total number of two-qubit gates needed in the robust 
CNOT construction is 

nf = lOWr 1 , i= 1,2,... , (14) 

again assuming the same N T for each level of concatena- 
tion. Additional re-isolating and concatenation increases 
the number of single and two-qubit operations required. 
We may be able to reduce the total number of single 
qubit operations by compounding gates however this is 
not possible for the two-qubit operations. The process of 
re-isolating the Ising component of the two-qubit oper- 
ation slices the rotation into many smaller rotations, as 
well as adding extra operations. The limit to how finely 
we can slice will be decided by the strength of the ex- 
change coupling. To perform small two-qubit rotations 



we may require a small J such that the operations evolve 
slowly enough to be within the realms of experimental 
pulse timing control. The viability of using multiple con- 
catenation for constructing robust two-qubit gates lies 
in tenuous balance between the ability to perform the 
large number of operations required quickly, and ade- 
quate pulse timing control over the small two-qubit ro- 
tations. The strength of the exchange coupling of our 
system will determine whether these conditions can be 
satisfied. 



B. Gate time 

Each level of concatenation increases the time taken 
for the robust CNOT operation significantly. In a work- 
ing quantum computer this may be problematic as the 
decoherence time of the system sets an upper limit on 
how long operations may take. For phosphorus donors 
in Si the coherence time, T2, of donor electron spins has 
been measured to be T2 > 60 ms at 7 K— . We calcu- 
late the total time taken for the robust CNOT gate for 
various implementation levels based on gate times using 
global control methods^. The results for this appear 
in Table HI As in Ref . [TH, we assume that single qubit 
rotations by an angle it take 40 ns to perform as does 
the Hadamard gate. We also assume that two-qubit ro- 
tations by 7r/4 take 1.96 ns if the coupling strength is 
given by Jq = 0.132 fieV, taken from the calculated un- 
biased exchange data*. Actual time will decrease under 
the application of a J-gate bias^, however, we assume 
a worse case scenario here. 



TABLE I: CNOT gate times for various pulse implementation 
levels in the electron spin solid-state quantum computing ar- 
chitecture. 



implementation 
level 


gate times (ns) 
single qubit two-qubit 


total 





180 


3.92 


183.92 


1 


716 


35.28 


751.28 


2 


53256.80 


2544.08 


55800.88 



As Table U demonstrates, operation time grows appre- 
ciably with concatenation. Furthermore, Fig [5] shows 
that the success of the robust CNOT gate is dependent 
on how accurately we can estimate the exchange coupling 
strength based on expectations of the fabrication process 
alone. In such an uncharacterised system we have shown 
that a sensible choice can be made based upon the tar- 
get separation, yielding Jo. Large variations in the ex- 
change interaction strength due to donor misplacement, 
and the additional time cost for multiple concatenation 
means composite rotations alone can not always guaran- 
tee a feasible, robust CNOT gate. However, we will now 
show that composite pulses at the lowest level coupled 
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with a systematic two-qubit interaction characterisation 
procedure allows for precise CNOT gate construction. 



III. THE ROLE OF TWO-QUBIT 
HAMILTONIAN CHARACTERISATION 

Using a combination of system indentification and 
composite rotations, we may construct a high fidelity 
robust CNOT gate. Whilst many methods of system 
identification exist, we choose the procedure of Hamil- 
tonian characterisation because it provides direct knowl- 
edge of the Hamiltonian (which we require) in an effi- 
cient manner. This approach strikes a balance between 
the need for multiple concatenation and precision Hamil- 
tonian characterisation, and may be particularly useful 
for systems whose Hamiltonian parameters require re- 
characterisation over time due to drift. 

Recent work shows how characterisation of a two-qubit 
Hamiltonian can be achieved via entanglement mapping 
of the squared concurrence relatio n 16 ! 17 . The identifica- 
tion of the Hamiltonian coefficients amounts to determin- 
ing the oscillation frequency of this entanglement func- 
tion for different input states. The only requirements are 
an accurately characterised Hadamard gate and measure- 
ment on both qubits. An important result from the work 
in Ref. [lil is the fractional uncertainty in a frequency 
determination 



> 



(15) 



where, N t is the number of discrete time points at which 
N e projective measurements are made. An equivalent re- 
sult can also be found in the earlier work of Huelga et al. 
in the context of Ramsey spectroscopy^. To accurately 
determine the frequency, the time over which the system 
is observed, i b, should be maximised, however this pro- 
cess is limited by the decoherence time of the system. 
An accurate frequency determination is still possible in 
the presence of decoherence by allowing t D b to be rela- 
tively large and performing two measurements at Nt time 
points. The uncertainty in the frequency can then be re- 
duced by evolving the system for a suitably long time 
before measuring at two final time points. This process 
is repeated N e times to estimate the phase of the oscilla- 
tion. The total number of required measurements is then 
TV = 2(N t + N e ). Characterising the system in this way 
results in the scaling of Eq. [T5] 

To characterise the Heisenberg Hamiltonian with 
isotropic couplings requires determining the oscillation 
frequency of three different input states, meaning N = 
6(Nt + N e ) total measurements are needed. The frac- 
tional uncertainty in the characterised exchange cou- 
pling, J c , as a function of N for a given Nt is 



5Jc 

Jc 



> 



Ay/6 



To illustrate the effect of composite rotations we consider 
a modest amount of characterisation by choosing Nt = 
10. Increasing the number of time points results in higher 
precision characterisation. 

In an uncharacterised system we assumed the coupling 
between donors, Jo, to be determined by the target donor 
separation. Donor misplacement as a result of fabrication 
uncertainties lead to variations in the coupling strength, 
J, from the target J . We have seen how the robust 
CNOT gate for an uncharacterised system performs in 
Fig. [SJ We now consider the performance of a robust 
CNOT gate in a characterised system. 

Characterisation of the Hamiltonian can be performed 
to any level of precision at the expense of extra mea- 
surements, with the uncertainty given by Eq. 1161 In a 
characterised system, the estimated coupling strength is 
set to the characterised coupling strength, J c (with un- 
certainty bounds ±<5J C ), rather than Jo. The fractional 
error in this case is 



Ac 



J 

T 



i. 



(17) 



f ~ NWN-6N t ' 



(16) 



where in general the characterised coupling strength, J c , 
will be much closer to the true value of J, than the tar- 
get value, Jo is to J. This means higher fidelity can be 
achieved using fewer levels of concatenation. 

Given that the total gate time increases so sharply with 
increased concatenation, operating with a single appli- 
cation of the composite rotations is preferential. For a 
one site deviation from the target separation, we show 
the resulting CNOT fidelity as a function of pulse im- 
plementation in a system characterised to the 10% level 
(5J C /J C = 0.1) in Fig. O Characterisation to this level 
would require at least 156 measurements assuming the 
previous parameters. We take J c f» 0.9J to be the char- 
acterised value of the exchange coupling strength, as it 
corresponds to an extremal bound value. The results 
in Fig. [6] demonstrate that it is possible to construct a 
very high fidelity CNOT gate using one level of robust 
pulsing, provided a suitable amount of characterisation 
is first performed. 

The total number of characterisation measurements 
needed to achieve a given fidelity can also be determined 
as a function of the implementation level. These results 
appear in Fig. [7] In reality the fidelity may be substan- 
tially higher than the results of Fig. [7] indicate, as they 
provide a lower bound for the corresponding number of 
measurements. These results show the clear benefit in 
using a single level of composite rotations and charac- 
terisation to construct a robust CNOT gate. The im- 
provements expected beyond this do not seem to warrant 
concatenation. 

Any quantum computation proposal requires that 
many operations be performed within the dephasing 
time, T2, of the system. The 10 -4 level is widely assumed 
to be the fault-tolerant threshold for both environmen- 
tally induced and systematic errors^, however more rig- 
orous bounds^ recently calculated, suggest it could be 
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FIG. 6: Exchange interaction strength as a function of donor 
separation along the [100] direction, showing a large variation 
in the coupling strength with donor misplacement (dots indi- 
cate actual site separations). For an uncharacterised system 
the coupling is set to the fabrication target Jo, with the ac- 
tual placement giving coupling J. The resulting CNOT error, 
for a one site deviation (Ao ~ —0.49) from the target 
separation can be seen on the inset plot as a function of im- 
plementation level. In the characterised system the coupling 
is set to J c . The CNOT error for a system characterised to 
the 10% level (8J C /J C = 0.1), taking J c « 0.9J (A c = 0.1), is 
shown as a function of implementation level inset also. Note 
that all curves are included purely to guide the eye. 



closer to 10~ 5 . Figure [7] shows that it is possible to con- 
struct a CNOT gate to this precision level in the presence 
of significant fabrication induced uncertainties, using ei- 
ther multiple concatenation of the composite rotations or 
a combination of the composite rotations and a modest 
level of characterisation. 

Assuming the system has been characterised to a mod- 
est level beforehand, we now show that in order to re- 
main below the threshold for environmentally induced 
errors also, the robust CNOT should be constructed us- 
ing a single application of composite rotations and char- 
acterisation. In Fig. [5J these results are shown for a 
system with an unbiased J-gate, J(V = 0), based on 
the 60 ms dephasing time in isotopically pure 28 Si at 
7 K, and for characterisation to the 10% level, again as- 
suming the extremal bound value of J c ~ 0.9J. In a 
biased system, the exchange coupling is stronger. Cal- 
culations suggest that for donors separated by ~ 20 nm 
in the [100] direction, a 1 V bias applied to the control 
gates can strengthen the coupling by over two orders of 
mag nitudo 4 .^. A robust CNOT gate comprising charac- 
terisation as described above could therefore operate at 
close to the 10 -7 level for environmentally induced errors. 
Performing additional measurements to characterise the 
system to the 1% level would lower the systematic error 



FIG. 7: CNOT error, (1 - f), as a function of the total 
number of characterisation measurements required to achieve 
a given fidelity for various implementation levels. The results 
demonstrate the usefulness of combining composite rotations 
with Hamiltonian characterisation when constructing a robust 



CNOT gate, 
shown. 



Threshold reference line at 10 error rate is 



level to well below 10~ 7 also, bringing it well within more 
rigorous threshold bounds^. 

For systems whose Hamiltonian parameters are not 
well known due to fabrication uncertainties, or may drift 
over time, this is an important result, suggesting that op- 
crating the CNOT gate in this way can guarantee that the 
error rate remains below the fault-tolerant error thresh- 
old. For the case of Si:P quantum computer architec- 
tures Fig. [8] suggests that this may be fabrication un- 
certainties within up to six sites of the target site, or 
~ 6.5 nm in the unbiased case, however in the J-gate 
biased case this allowance may be much greater. The 
trade-off for operating in this manner is the need for pe- 
riodic rc-charactcrisation, however the cost of this should 
be minimal as the number of required measurements is 
small. 



IV. CONCLUSIONS 

The performance of a robust CNOT gate constructed 
using two-qubit composite rotations has been examined. 
Multiple concatenation of the composite rotations re- 
sults in a high fidelity CNOT gate provided the frac- 
tional uncertainty in J lies within the correctable range. 
Large variations in the exchange interaction coupling 
with donor separation means this is not always the case. 
Furthermore, multiple concatenation of composite rota- 
tions requires long overall gate times with respect to the 
dechorence time of the system and results in gate oper- 
ation which exceeds the current error threshold required 
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^ 10" 
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10" 5 



10" 



uncorrected 



level comp. 



error threshold rcf. 



1 st level comp./10% char. 



Misplacement by: 
1 site O 
2 sites x 
6 sites □ 
9 sites + 



level comp. 



iL 



1Q -b 1Q -i, 1Q -4 

7cnot/T2 

FIG. 8: CNOT error, (1 - T), as a function of the total gate 
time for an unbiased, J(V = 0), system (T2 = 60 ms as- 
sumed). Results are shown for a range of separations in the 
[100] direction, larger than the targeted 20.634 nm separation. 
We consider various CNOT gate constructions, namely an un- 
corrected CNOT, one constructed from both a single and two 
applications of composite rotations and finally a CNOT con- 
structed using composite rotations in conjunction with char- 
acterisation to the 10% level (<5J C /J C = 0.1) taking J c ~ 0.9J. 
Only for this final method have more than two sites been 
included as for other methods results will clearly be worse. 
Threshold reference lines at 10 -4 error rates are shown. 



for fault-tolerant quantum computation. As an effective 
fix to this problem, we demonstrated how, in a system 
with large variations in the qubit coupling strength, a 
high fidelity CNOT gate which operates below this error 
threshold can be constructed from a single level of com- 
posite rotations in conjunction with Hamiltonian charac- 
terisation. 
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